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The effect of electric fields on lipid membranes 
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Contrary to existing theoretical models, experimental evidence points out that electroporation 
(membrane defect formation under external electric fields) starts to occur within the range of trans- 
membrane voltages that cells may routinely experience, curiously, just above the range of trans- 
membrane voltages involved in neural signal transmission. Understanding the underlying principles 
of electric fields-lipid membrane interactions seems to carry a great biological importance. 

An argument is presented toward understanding the theoretical aspects of electroporation by 
using the DLVO theory, which has not been recognized previously in the context of electroporation. 
Further, the dispersion interactions (with its quantum nature), of the double layer counterions and 
membrane lipid molecules over the Stern layer are emphasized. The sign of these forces is such that 
they compress the membrane. A parallel is drawn to the theory of thin films. The argument is that 
the external electric field breaks the symmetry of the disjoining pressures on both sides of a lipid 
membrane, resulting in a protrusion of only few lipid molecules. That compromises the membrane 
stability on a nanoscale and makes it traversable to ions. The presented estimate based on these 
arguments is consistent to electroporation experiments and existing numerical simulations. 
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I. INTRODUCTION 



Foam films and membranes have been studied for over 
100 years. They are a bilayer formation with or without 
an aqueous core respectively, in which the hydrophilic 
groups are facing the aqueous solution. In the case of 
membranes the hydrophobic ends are in opposite direc- 
tion facing each other, and in the case of foam films, 
soap bubbles, for example, the hydrophobic ends are ex- 
posed to air. In both, foam films (soap lamellas) and 
double layer membranes (lipid membranes), an equilib- 
rium on the order of molecular dimension thickness (nm) 
can be reached, creating black films and membranes [l|. 
The biomembranes play a vital role in the activities of 
organelles and cells, and they are very closely related 
to this group of bilayer formations. Their main function 
varies, but one of the essential roles is to provide an excel- 
lent barrier, allow a very complex and selective molecular 
signaling and transport, including propagation of neural 
signals. 

Beside diverse composition and thickness of the bilayer 
and its exposition to air or aqueous medium, they all 
exhibit more or less elastic stability. Under mechanical 
or electrical stress they rupture (lyse) and in some 
occasions they may exhibit reversibility of the defects 
(reversible breakdown). The reversible breakdown has 
been observed in the case of soap films exposed to elec- 
tric fields for example 0] and is a more common feature 
of lipid membrane electroporation. This indicates an 
underlying mechanism how these disperse systems can 
be analyzed and described. 
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A. Electroporation 

The electrical breakdown and reversible breakdown of 
biological membranes has been subject of intensive exper- 
imental and theoretical studies. However, the effects the 
electric fields have on membranes, and the actual mech- 
anism of electrically induced instabilities is still not sat- 
isfactory understood. Some models are concerned with 
more macroscopic reasons for instability as lipid mem- 
brane thinning, related to electrostriction , undulation 
and elastic properties of the membranes [!, H, [H, [f| . Other 
models, as electroporation models, find the explanation 
in formation and expansion of tiny pores. This process 
of electrically creating aqueous pores in the membrane, 
large enough to be traversable to ions and molecules is 
termed electroporation. The effect is believed to be in- 
trinsic to artificial planar bilayer membranes, as well as 
to any other biological membrane of a cell or an organelle. 
Thus, the electroporation process may strongly increases 
the ionic and molecular transport across a biologically 
relevant membranes, redistribute membrane's lipid con- 
tent and cause many biologically significant effects. 

The process of electroporation has been routinely used 
in cell biology to load cells with small molecules (drugs), 
and large molecules (as proteins and DNA) [7[. 

Experimental evidence over many decades Q sug- 
gested that the electroporation process is a stochastic 
process, with pore size evolution and occurrence of re- 
versible or irreversible breakdown [ij [l(| ■ It is believed 
that the pores completely dominate electrical conduc- 
tance, but they occupy very small fractional area. 

Bilipid membrane thickness is much smaller than the 
wavelength of light so the process of electroporation is 
most commonly observed through change of electrical 
conductance of the membrane (believed as a result of 
electrically induced pores). Recent accurate measure- 



2 



merits have recorded a change in conductance due 
to a single defect (pore) at transmembrane voltages as 
little as w 250 mV. The pore size estimates based on the 
change in conductance indicate a nanometer size defects, 
and a lifetime in the range of milliseconds. 

The theoretical models of electroporation are roughly 
concerned with two distinct aspects, the theory of pore 
creation and the theory of pore evolution. 

The pore evolution concepts have roots (l2l . fl3L IbH ] in 
the work of J.B. Zeldovich, from 1943. Based on exper- 
iments and theoretical analysis, the Chizmadzhev group 
in 1979, has described the initial electroporation theory 
for artificial planar bilayer membranes in a set of seven 
papers Q . The use of Smoluchowski equation in the con- 
text of pore size evolution was given by d, . 

While the pore evolution and expansion models seem 
to be well established, the initial steps of pore formation 
are less certain. A central part of almost all electro- 
poration models is the critical permealization threshold 
needed for pore creation. A phcnomcnological, Arrhcnius 
like model of transmembrane voltage de pen dent energy 
barrier was independently proposed by |l6l. More 
recent electroporation simulations [HI, Il9l] use different 
pore evolution concepts, but a common Arrhenius en- 
ergy barrier of w 45 — 50 kT for pore creation [20| , which 
is in agreement with the experimental values of critical 
voltages causing electroporation. The existing theoretical 
models, related to membrane elasticity, seem to predict 
critical voltages of several Volts causing membrane in- 
stability. That exceeds the experimental lower threshold 
values (w 250 mV - 1 V) by far. Model modifications 
have been proposed (2lj to avoid this discrepancy, but 
it seems that mechanism of membrane's electrical break- 
down is still not satisfactory understood. 

Besides all this advances, it seems that many questions 
related to the physical mechanisms of the electroporation 
are still not well understood [22j , and more basic research 
is needed. 

Based on the accumulated experimental observation 
and the recent measurements, specifically, [Tlj ]. and 
molecular dynamics simulations [23J], a few main elec- 
troporation aspects, relevant for this paper, can be 
summarized. 

- Relatively low transmembrane voltages of « 200 mV 
can cause electroporation. 

- Poration is a local effect (just few lipids involved). 

- Pores are nanometer size defects. 

- It takes a short time (ns) to create a pore. 

- It takes a long time (ms) to destroy a pore. 

- The process is sensitive to concentration and pH. 

- The process causes electroporation asymmetry and 
lipid redistribution on the cell membrane (flip-flop). 

This paper, further mainly addresses the theoretical 
aspects of pore creation. The other aspects of the elec- 
troporation will be addressed in a future work. 

To put all this into context, a few basic aspects of lipid 



membranes and mean field theories of ion interactions 
with membranes, are briefly noted. 



B. General about lipid membranes 

1. Lipid membrane fluctuations 

Based on spectroscopic data a time scale for reorien- 
tation of chemical bonds, molecular segments, molecules, 
or assemblies of molecules can be placed. In general, the 
correlation time increases with increasing size of moving 
segments. The time scale of lipid motion in biomem- 
branes and their approximate correlation times cover al- 
most 15 orders of magnitude [24J. Protrusion, the head 
groups sticking in or out of the membrane plane is on 
the order of 10~ 9 s. Rotational diffusion of the head 
group or wobble is on a scale of 10 -8 s. The lateral dif- 
fusion (in the plane of the membrane) is on the order of 
10~7 s, with a diffusion constant of D w 10~ 12 to 10~ n 
m 2 /s, and lipids spanning distances of w lO^im in 1 s. 
The large scale collective deformations (curvatures) from 
the plane of the membrane - undulations are on a time 
scale of 10 -6 s to 1 s. They are thermally driven, and 
membrane's waves like motions correspond to very large 
wavelengths. Finally the jumps of lipids from one side 
on the membrane to the other or the so called flip-fop is 
dependent on molecule's polarity, and is on the order of 
10 3 s to 10 4 s. 



2. The pressure profile in lipid membranes 

Protrusion of the lipid molecules is closely related to 
the lateral pressure profile. Most of the biologically rele- 
vant lipid bilayer membranes are made of lipid molecules 
with two hydrocarbon chains [25|. The lateral pressure 
profile of these lipid bilayer membranes is symmetric with 
respect to the middle of the membrane. Going towards 
the middle of the bilayer [25j , there is initially headgroup 
repulsion, and then interfacial attraction. Coming to- 
wards the two hydrocarbon chains, repulsion is encoun- 
tered again, since hydrophobic hydrocarbon tails dislike 
each other's presence. All of these pressures are summed 
up to zero. 

The value of these pressures is hig h. Various molecu- 
lar dynamics simulations [26|, [2?], l28l l29j indicate values 
in the range of hundreds of atmospheres. A simulation 
estimate using lattice statistical thermodynamic calcula- 
tions of lateral pressure profiles [30( give a value of over 
400atm. A lower threshold value of « 400atm will be 
used throughout this paper as a approximate threshold 
of membrane instability. 

Molecular dynamics simulations of electroporation [23| 
also indicate high values of pressures. The forces of tens 
of pN on the area of a lipid headgroup (« 0.6 nm 2 )[25|. 
result in fairly large pressures of few hundred atmo- 
spheres. 
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C. Mean Field Theory of Ions near lipid membrane 



II. DLVO THEORY AND MEMBRANES 



The surface charge density of lipid membranes mainly 
comes from the lipid head groups that are exposed to the 
aqueous solution. The electric field of this surface charge 
redistributes the ions in the electrolyte due to the local 
concentration (n) and potential (ip) gradients that can 
be described according to the familiar equation 



dp d 
dt dx 



<9p ze dip 
die + ~kT P ~dx~ 



(1) 



where planar symmetry is assumed, x being perpen- 
dicular to the membrane, and the charge density p — 
^2 Pi — ~Yli n i z i e i s a sum over all the concentrations of 

i i 

ions rii, of valence z%. At electrochemical equilibrium we 
have 



An important contribution to the problem of stability 
of thin films was given by the works that started 60 to 70 
years ago by Derjaguin, Landau, Verwey and Overbeek 
(DLVO theory). The behavior of the system of electrolyte 
column sandwiched between any two surfaces is a conse- 
quence of the interplay between repulsive and attractive 
(van der Waals) forces on the interface. The balance of 
these forces and the thickness of the film correspond to 
an interaction energy minimum. 

In 1939, Derjaguin originally formulated the concept 
of a disjoining pressure for thin liquid films and verified 
experimentally its existence [l3| Changes in the 

interfacial region that generate the disjoining pressure in 
a thin liquid film originate from intcrmolccular forces. It 
is customary to separate the various contributions of the 
disjoining pressure into different components, e.g., 



Pioe 



(2) 



where pio is the bulk charge density far away from 
the membrane (ip (x) ps 0). Then the potential of this 
ion distribution at equilibrium is given by the Poisson- 
Boltzmann equation 



dx* eeo 2^ z ^° e 



(3) 



Rearranging equation ([3]) it can be seen that the fol- 
lowing expression is constant with x. 



d 

dx 



E 



n kT — — ( — 

nihL 2 [dx 



= -[P(,)]=0. (4) 



The expression in the brackets (pressure) shows that 
the entropy of mixing (the osmotic pressure) of the elec- 
trolyte is influenced by the electric field from the mem- 
brane surface charges. As a result, entropically unfavor- 
able, counterions are more densely distributed against 
the membrane, but in such a way that the total pres- 
sure is constant in any transverse cross section from 
the membrane. The expression Q is dependent on the 
ion concentration, and by applying an additional electric 
field, there will be ion redistribution. It should be men- 
tioned jTsj] that it is not necessary to resort to Poisson- 
Boltzmann equation and any assumptions leading to it, 
to obtain expression In most general case we would 
have 



d 

dx 



p(ip)dijj 



(to 



dip 
dx 



dx 



[P(x)] = 0. (5) 



n = n + n w + n s 



n. 



adh 



(0) 



where the first two, the electrostatic (LI C ) and the 
dispersion (II W ) component, usually refer to DLVO 
theory. The other components as steric, short range 
structural forces, and adhesion forces play important 
role at molecular distances. The key assumption of 
equation Q is that the various contributions to the dis- 
joining pressure are additive. However, it is not always 
clear that this assumption is valid and it may lead, in 
some cases, to anomalous results. In general, thin film 
stability is influenced most importantly, by the electro- 
static (LI e ) and the dispersion (LT W ) components [l3| 

Further the disjoining pressure components of a single 
lipid membrane in solution are analyzed, by considering 
the case when one of the films is moved far away. 



A. The electrostatic component of disjoining 
pressure n c 

DLVO describes the interaction pressure of thin films. 
There the diffuse layers of the two charged surfaces over- 
lap, and cause interaction pressure between them. The 
hydrostatic pressure at some distance x from the mem- 
brane is given by the constant term in equation ([5]). The 
difference of the hydrostatic pressures in the bulk (ab- 
sence of the other surface) and in the interlayer is the 
electrostatic component of the disjoining pressure II e [l3| 



n„ 



Pd(x) - Poo(x). 



(7) 



Here Pd ix) is the pressure at point x in the interlayer 
when the two surfaces are 2D apart, and Poa{x) is the 
pressure at the same point x when the two surfaces are 
separated. 
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If the overlapping diffuse layers are identical, there is 
a symmetry plane at x — D where ip — ip m and E = 0. 
Considering the anions and cations, in a 1 : 1 electrolyte 
such as NaCl or KCl, and a Gouy-Chapman p(ip), the 
disjoining pressure from equation ([7]) will have a cosine 
hyperbolic dependence on the midplane voltage ipm, 



n e (ipn 



p(ip)dip = —2nokT 



.(8) 



For a single film (membrane), when the other film is 
far away, but D is still finite, the fact that the mid plane 
potential ip m is very small (tp m ~ 0) is used, and the 
following approximate expression (as derived by Verwey 
and Overbeek in 1948 [3l|) can be used. 



n e (tp m ) ~ -n 



(zey 
kT 



(9) 



The pressure is directed towards the membrane. 

Further, in this paper, arguments are presented that 
the expression © should be used to describe the inter- 
action of external electric fields with a lipid membrane 
in electroporation models and theory. The argument for 
this follows the analogy that the presence of electrode's 
electric fields has a similar effect as another charged sur- 
face being present far away. 



B. The dispersion component of disjoining 
pressure II W 

The van der Waals dispersion forces are caused by po- 
larization due to fluctuating electromagnetic fields. Be- 
sides being dependent on geometry, they are also fre- 
quency dependent, and they relate to the refractive index 
and dispersion. 

It can be shown, for example (32l . |33| that van der 
Waals interaction of two neutral polarizable particles 
can be attributed to correlations in the fluctuating 
vacuum electromagnetic field, that induces fluctuating 
dipole moments in the particles, and the correlation of 
these moments is the der Waals interaction. The van 
der Waals interactions can be derived from the Maxwell 
equations and the additional (quantum) assumption that 
each mode of the electromagnetic field has a zero-point 
energy \l 



\hoj. 



Calculating the dispersion forces, even for simple 
geometry, is not an easy task and few aspects like 
nonadditivity, relativistic retardation and the sign need 
to be addressed. 

In rarified media like a gas or a solute in a solution, 
non-pair wise interactions are small and additivity holds, 
but in condensed media, the dispersion forces are not 



additive, and in general, the force between two molecules 
depends on the presence of other molecules. These 
non-pair wise interaction energies may be attractive or 
repulsive, depending on the geometrical arrangement 
of the particles involved, and are important to be 
considered if the intermolecular spacing is small. The 
Lifshitz formula and theory [12, HH, [32j used here to 
calculate the Maxwell stress tensor is macroscopic, and 
avoids the problems associated with non-additivity of 
the dispersion forces. 

When the intermolecular spacing is large, the rel- 
ativistic retardation, due to the finite speed of light, 
needs to be included. This aspect plays an important 
role in colloidal systems and reduces the strength of 
van der Waals interaction [34[. In the limiting case of 
perfect conductor, including retardation, the Lifshitz 
formula reduces to the Casimir force, which is a macro- 
scopic manifestation of the van der Waals force. In 
either of these cases van der Waals interactions may be 
regarded as a consequence of the fluctuating vacuum 
electromagnetic field |32l . [33j . For the distance between 
the interacting objects considered in this paper, the 
relativistic retardation can be ignored. 

Regarding the sign of the dispersion forces, they can be 
positive or negative depending on the situation and the 
geometry. Form the de Boer-Hamaker theorem [l3|, H|| , 
it follows that the van der Waals forces between two par- 
ticles of fixed shape and of the same material embedded 
in a fluid is always attractive. If particles have different 
composition, the resulting force can be repulsion. This 
is because in the interaction between the two bodies, an 
interaction of fluid bodies of same size has to be included. 

Considering all this, an expression and the sign for 
the dispersion component of the disjoining pressure II W 
in the context of electroporation, is further presented. 
This expression is given for a solute particles and a 
substrate, both in medium. The assumption is that the 
fluctuating electric field in a dilute solution or gas, acts 
on each dilute species without distortion of the field 
by other solute or gas molecules. For such media the 
Lorentz-Lorenz corrections are unimportant, the Lifshitz 
formula can be simplified, and e of the dilute-solution 
(dilute-gas) limit has the dielectric response strictly 
proportional to the number of solute (gas) molecules 

Hall! S3- 

The interaction energy of one particle (e s ) of radius a, 
at a distance d away from a substrate (ehc) over a medium 
(e w ) is given as [36| 



kT ^2. 



2d 3 ^ 47T£ V 
n— 



6} 1C 6 V 
^hc + C v 



(10) 



where /3/47re w = a 3 (e s — e w )/(e B — e w ) and all three 
dielectric permittivities e(i£) are dependent and summed 
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over all of the imaginary frequencies i£„. The prime over 
the summation indicates that the zero frequency term 
(n = 0) is multiplied by 1/2. The force per unit area, 
exerted by each ion on the substrate would be a deriva- 
tive of equation (|10p . Then the dispersion component of 
the disjoining pressure for the number of particles in a 
segment Ax « a of unit area, parallel to the substrate, 
given in terms of the particle's concentration n(d) is 



n w (d) 



Axa 3 3h f (e s — e 



(Che 



d 4 2 J (e s + 2e w ) (e hc + e«) 
o 



dc (ii) 



For the dispersion part of equation (JTUJ) , due to the 
fact that the electronic polarizabilities are dominant, the 
typical absorption frequencies £ are close to the UV re- 
gion. Then the pre factor kT and the sum over discrete 
frequencies in equation IjlOp , are replaced by h and an in- 
tegral over them [25| , thus revealing the quantum nature 
of equation (|11[) and the role played by the zero-point en- 
ergy. Calculating the integral in equation (|11[) . requires 
the functions e(i£), to be known. Unfortunately, com- 
plete data of the absorption spectra of most materials 
are not available. An estimate of the integral is given in 
the Appendix. 

From equation {TT]) it can be seen that the sign of the 
pressure II W is determined by the dielectric permittivi- 
ties for given frequency. The particles will be attracted 
towards the substrate if e s > e w and £hc > tw- 



ill. ANALYTICAL AND NUMERICAL 
ESTIMATES OF THE DISJOINING PRESSURES 

An 

A. Analytic estimates of the equilibrium potential 

From equation Q, noting that at the surface E s = 
u/eeo and in the bulk p + = p_, we get Grahame's equa- 
tion that relates the known surface charge density a with 
the unknown surface potential ip s in terms of inverse hy- 
perbolic functions. 



/ 2kT J 
ip s = ash 



(12) 



From it and from equation ^ , the surface charge den- 
sity p s can be calculated. In the case of thin films, the 
mid-plane (where E = 0) is very useful, but this is not 
applicable in the case of a single membrane. 

The potential at any distance x from the surface in 
Gouy-Chapman Theory (high potential) [25[ is given as 



ip(x) 



AkT 



(13) 



where 7 = th(eijj s /4kT). The Debye length 1/re ap- 
pears as the characteristic decay length of the potential, 
same as in Debye-Hiickel theory. 



1. Numerical simulation 

A coupled simulation of equations (JTJ) and © over 
time of 1 ps and simulation space of 0.25/im was done. 
The bulk ion concentration was considered to be lOOmM 
which corresponds to Debye length of ps lnm. Micro- 
scopic details on the membrane were not considered. It 
was taken to be a 5nm thick dielectric slab with dielec- 
tric constant of e = 2. On both sides of the membrane 
a Stern layer of thickness 1.8nm, approximately the size 
of a water molecule [25[ , not accessible to diffusing ions, 
was considered. The dielectric constant in the Stern layer 
was taken to be the same as the one in the water, e = 80. 
The analytical expressions (TT2"]) and (TIB"]) were used in the 
simulation of equations U) and ((3|) to find ip (x) and from 
it, the initial ion distribution p (x) near the membrane. 
This initial ion distribution was used later in the case 
when an external electric field was applied. 

The simulation gives a nice agreement and compara- 
ble results to Gouy-Chapman theory (high potentials), 
and the Debye-Hiickel theory (low potentials). Under an 
external electric filed there is an ion concentration stabi- 
lization near the membrane within the membrane charg- 
ing time of w 1/iS. This also corresponds to the standard 
membrane specific capacitance of c m = luFcnU 2 . The 
analytical estimate of the ion concentration change based 
on the specific capacitance equation (|16[) gives a value of 
An s» 14.4mM for each increase in the applied voltage of 
A^ m = 50mV. For the run time of the simulation (lfis) 
the ion concentration change was An ps llmM. 

The type of ions, the bulk ion concentration, and the 
applied voltages used in the simulation were chosen to 
resemble the experimental situation of [Tlj . Besides the 
agreement, the data from the simulation was used just to 
illustrate the voltage profile near the membrane in figure 
@. In all other estimates, standard analytical values 
were used. 



B. Estimate of the gradient of the disjoining 
pressures 

1. The gradient of the electrostatic component of disjoining 
pressure Ail c 

The electrostatic component of this pressure gradient 
over the membrane can be calculated using equation J9J , 
where ?/; m is the value of the potential far away from 
the membrane. The applied electric field manifests it- 
self by establishing Aip m across the membrane without 
substantially affecting the side potentials [37( near the 
membrane. The side potentials ip s depend heavily on 
membranes surface charge and ip S 2 = tp s i+Ai() m as shown 
in figure ([2]). Then relative to the left side potential, the 
difference in the electrostatic component of the disjoining 
pressure on both sides of the membrane is 
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An c (v> m ) = n loft - n right = n ^-Av£, (w) 

where A?/>m is the transmembrane voltage. This ex- 
pression has not been used yet in the electroporation 
models and theory. Besides its voltage dependence, away 
from the membrane, it also reveals its bulk ion concen- 
tration dependence that is often missing in the electro- 
poration models. 

In electroporation experiments [ll[ , it has been noticed 
that relatively small transmembrane voltage of few hun- 
dred mV cause electroporation, in contrast to membrane 
elasticity model estimates of few volts. The expression 
(H3J) indicates that Aip m « 300 mV gives AII C over the 
membrane instability threshold of ~ 400 atm. The cal- 
culated values of An e are shown in figure Jl}. 



d<Ax<d + a. (17) 

According to [251 ] . the hydrated ion radii are a-^ a + = 
3.6A, and ok+ = 3.3A. Finally, taking the lower limit 
of equation (fT7|) , for K + ion concentrations of An = 5 x 
14.4 mM, corresponding to ?/> m = 250 mV we will have 
a value for the lower limit of AII W = 113 atm (146 atm 
for Na + ions) . The values for II W (for ?/> m = 250 mV) are 
comparable to n o , and their sum is in the range of the 
membrane instability threshold of « 400 atm. Estimates 
for wider range of transmembrane voltages are shown in 
figure (U). 

From equation (|15[) it can be seen that the contribution 
to II W of ions at larger distances diminishes rapidly. 

C. General Electroporation Conclusions 



2. The gradient of the dispersion component of disjoining 
pressure AI1 W 

Similar to AII C , due to the externally applied field, 
the asymmetry in the number of ions on both sides of 
the membrane (ni c f t = n oq + An and n r i g ht = "cq 
\ ; will also cause a pressure gradient in the dispersion 
component II W . As indicated from equation (fTTj) van der 
Waals dispersion force is proportional to the hydrated 
ion concentration n{d), in a segment Ax of unit area, 
parallel to the membrane. The dispersion force between 
the ions and the membrane hydrocarbons will give rise 
to pressure gradient AII W due to the difference in the 
number of ions on both sides of the membrane. 



AIL 



2An a 



, Ax 
~dF 



(15) 



where / (in units of energy) is the integral over the 
imaginary frequencies in equation (|11J) . An estimate of / 
is given in the Appendix. Further, AII W can be estimated 
for different types of ions. 

An is proportional with the accumulated charge and 
an estimate of An that goes into equation (|15[) can be 
made through membrane's specific capacitance c m , 



An = 



1 A^ n 



2edx 



■(l-e-*/^), 



(16) 



and for times longer than membrane's charging time 
t T c h all of the An accumulated ions will be next to 
the membrane at a distance of the Stern layer from it. 
The above expression gives An as 14.4mM w 0.2% on 
one side, for every Aip m = 50mV. When the An ions 
are being distributed next to the Stern layer (d — 1.8 A), 
which is about the size of the water molecule diameter 
[25l | , Ax is in the range of values 



Based on the DLVO approach, the main electropora- 
tion aspects now can be interpreted, and few conclusions 
can be drawn. 

It is been experimentally noticed [ll| that relatively 
small values (~ 250mV) of the transmembrane voltage 
cause electroporation, which has not been satisfactory 
explained yet. Arguments using standard DLVO theory 
(equations (|14p and (|15p) were presented here, that indi- 
cate how a relatively low transmembrane voltage results 
in large disjoining pressures (hundreds of atmospheres) 
on the membrane eventually causing electroporation. 

Based on membrane's electrical conductivity during 
electroporation, experimental evidence strongly suggests 
the existence of nanometer size defects. This makes elec- 
troporation a local effect that involves rearrangement of 
only few lipid headgroups. The timescale for creation of 
these defects is of the order of nanoseconds. The only 
lipid membrane fluctuations of this spatial and temporal 
scale are the protrusions. This suggests that the dis- 
joining pressure is affecting the protrusion of the lipid 
headgroups, and is working against their lateral pressure 
profile. This also suggests that the external electric field 
compromising the membrane's stability can be localized 
too. 

It takes a short time to create a pore. In contrast it 
takes a long time for pore destruction (pore lifetime is 
rs ms [11]). DLVO gives an insight why this might be 
the case. Once a nanometer size pore is filed with wa- 
ter, it resembles more the situation of thin films with 
overlapping interfacial regions. Van der Waals forces re- 
verse sign for this situation and prevent the closing of 
the pore on nanosecond timescale due to double layer and 
hydration repulsion, which in the case of thin films, dom- 
inate phospholipid membranes as they approach within 
0.1 to 0.3nm. The interfacial regions of these double lay- 
ers would depend on the geometrical shape of the pore, 
which at present is not known and is usually modeled as 
cylindrical or toroidal. Ion concentration inside the pore 
will also be different from the bulk concentration. 
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Analytical Estimates of the Disjoining Pressure for lOOmM KC1 (NaCl) solution 




Na + -> 
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FIG. 1: Analytical estimates of the electrostatic and dispersion components of the disjoining pressure. Equation (|14|l and 
(|15l) were used. The increase of the number of ions with the transmembrane voltage is at a rate of 14.4 mM per 50 mV 
(An = 14.4(Ai/> m /50mV)mM). The value of II ~ 400 atm is indicated as an approximate threshold for membrane instability. 
If instead of KCl a NaCl salt is used, II will be higher as indicated. Note that II W > II C for lower values of A?/> m . 



The DLVO expressions of the disjoining pressures are 
bulk concentration and ion type dependent. Expression 
© has not been used yet in the electroporation models 
and theory. Besides its voltage dependence, it also re- 
veals its bulk ion concentration dependence that is often 
missing in the electroporation models. 

The cell contains mostly K + while the extracellular 
space contains mostly Na + ions. If different ions are 
present on each side of the membrane, as in the example 
of a cell, the electroporation models predict that both 
the pore density and the transmembrane potential will 
be symmetric. But asymmetrical permeabilization of cell 
membranes has been experimentally observed [38, 39]. In 
light of DLVO, the difference in An to t may account for 
this asymmetry. The disjoining pressures on both sides 
of the cell will be different, since AII W will be different 
for K + and Na + counterions. The fact that the cell 
side facing the anode is being most permeable (38j is in 
agreement with the II W being larger for Na + than for 
K + counterions as seen in figure (fTj). 

Pores live long enough (ms) for lateral diffusion (i=s 
10~ 8 m in ms) to cause redistribution of lipids from one 
side to the other. 

Applying DLVO to present electroporation simulation 



models suggests slight modifications that need to be 
made. The phenomenological voltage dependent pore 
creation rate with the constant B that was crudely es- 
timated from Melinkov's experiments [2(j, needs to be 
related to the prefactor in equation (fl4|) and equations 
(|T5|) appropriately scaled by a sutable volume. This will 
make the pore creation rate dependent on the bulk con- 
centration and other material properties of the system 
contained in the constants, which so far was not the case. 
This will also slightly change figure (2) in [20( for the 
lower transmembrane voltage values, and relate the lag 
time with AIT W and the alignment of ions against the 
membrane. Since for Il w to reach its maximal value, An 
really needs to arrive within the Stern layer, next to the 
membrane, while Aip m and AI1 C do not change substan- 
tially. 

Some of these conclusions should be elaborated in more 
details in a future work. 



IV. CONCLUSIONS 

The standard DLVO theory, that has been used to de- 
scribe the stability of thin films, was applied to stability 
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FIG. 2: Profile of potential ip(x) across a membrane for different applied fields causing transmembrane voltage A^m of 
fOO, 200, 300 (indicated) and 400 mV. On both sides near the membrane, there is an additional jump of ip(x) ~ 100 mV in 
the Stern layer [25| ]. The Stern potential is not plotted on the left to indicate that A?/> m (and the applied field) have virtually 
no effect on the side potentials tp B Then it can be noted that ip E 2 = ip e i + Aipm. This fact is used in calculating AII C 

using equation (|14|l . This plot was obtained by simulating equations ([TJ and over time of 1 fis. Microscopic details in the 
membrane were ignored and it was considered to be a dielectric slab with surface charge density of a = 0.2 Cm~ 2 . 



of lipid membranes affected by electric fields. Electro- 
poration is caused by relatively small number of ions 
(» 1 — 2%), compared to the number of ions already 
present near the membrane. When external electric field 
is applied, it breaks the symmetry of the ion distribu- 
tion in the double layer on both sides of the membrane. 
This gives a rise to a pressure difference ATI in the dis- 
joining pressures on both sides on the membrane. The 
pressure gradient is causing protrusion of few lipid head- 
group compromising the stability of the bilayer. This idea 
was supported, by estimating the changes in the disjoin- 
ing pressure An to t, and comparing it with lower limit 
estimates of the bilayer's lateral pressure calculated from 
statistical dynamics calculations. In comparison with the 
existing electroporation models, the pore formation en- 
ergy per lipid volume used in these models corresponds 
to the Derjaguin's disjoining pressures. With respect 
to this, slight corrections to the existing electroporation 
models were also proposed. 

Finally, the dispersion interactions of the double layer 
counterions and membrane lipid molecules over the Stern 
layer were introduced in the context of electroporation 



and membrane stability. They seem to be dominant in 
the lower range of transmembrane voltages. They are ion 
size specific and they complement the electrostatic com- 
ponent of the disjoining pressure. Their effect on metal 
surfaces, for example, can be quite large, since expression 
(|A.5|) for metals is at least an order of magnitude larger. 
This can have significant effects on metal solution in- 
terfaces like electrodes in batteries, fuel cell membranes, 
electrochemical dealloying and similar situations. 

Experimental evidence points out that electroporation 
starts to occur within the range of transmembrane volt- 
ages that cells are routinely experiencing. This is just 
above the range of voltages involved in neural signal 
transmission (less than 200mV). Pores are created on 
a nanosecond timescale, but they persist for millisec- 
onds. Thus, it seems that the undesirable electroporation 
placed a limit on the range of transmembrane voltages 
involved in neural signal transmission. 

Main limitations of the standard DTVO approach pre- 
sented here are the facts that Poisson-Boltzmann and Lif- 
shitz theories are (continuous) mean field theories. Also 
not taken into account are few other effects, as the steric 
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FIG. 3: Applied electric fields result in ion concentration change of An(x) against the membrane, causing the indicated 
transmembrane voltages Ai/> m . With time more and more ions are coming to the closest possible approach to the membrane, 
which is the Stern layer (« 1.8A) [2^|. The simulation was run for time of 1/js and it indicates an ion concentration change of 
~ 11 mM in the volume left (right) of the membrane for each 50 mV increase in the transmembrane voltage. An analytical 
estimate for much longer times, using equation (|16[) indicates an accumulation of An = 14.4 mM in the immediate neighborhood 
of the membrane. These values are used to calculate AII W , using equation (|15[) and they are shown in figure (UJ). 



effects (due to finite ion sizes), image forces, solvation 
(hydration) forces, and discreteness of surface charges, 
between other effects. This may be a subject of future 
work. 



once 



All the data in the following estimate is from refer- 
251 ] . Calculating the integral in equation (fTTj) . for 



the dielectric permittivity in the visible and the ultravi- 
olet range that are dominant in the given situation, the 
following estimate is obtained. 
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APPENDIX: ESTIMATE OF AII W 

To calculate exactly the integral in equation (fTTj) . the 
functions e(i£), need to be known. Complete set of data 
of the absorption spectra of most materials are not avail- 
able. Reliable results can be obtained by several main 
bands in the absorption spectra, or to a first approxima- 
tion, the visible and the near ultraviolet bands [13j |. 



£hc = 1 
e™, = 1 



1 



n 2 — 1 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



where, according to [25|, the refractive indices are 
ri w = 1.33, rihc = 1.41, and the ion's dielectric permit- 
tivity e s is routinely approximated as that of a metal. 
The plasma frequencies of the free electron gas are £ c — 
3.0 x 10 15 i 



and £ se = 5.0 x lO^s" 1 . Using lA~2l lA~3l 
and IA.4I in equation (fTTj) , and solving the integral for £ 
gives the following value 
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specific ion types and concentrations. It should be noted 
x that this is just an order of magnitude estimate, con- 

3ft f (e s — e w ) (the — e w ) -20 sidering the imprecision of the quantities going into this 

^ — ~2~ y (e s + 2e w ) (the + £ w ) ^ ~~ ^ X ^ J(A-5) estimate, and the mean held approach. 

o 

This value of / is further used to estimate AI1 W for 
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